Dark-Matter Electric and Magnetic Dipole Moments 
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We consider the consequences of a neutral dark-matter particle with a nonzero electric and/or 
magnetic dipole moment. Theoretical constraints, as well as constraints from direct searches, pre- 
cision tests of the standard model, the cosmic microwave background and matter power spectra, 
and cosmic gamma rays, are included. We find that a relatively light particle with mass between 
an MeV and a few GeV and an electric or magnetic dipole as large as ~ 3 x 10 _16 e cm (roughly 
1.6 x 1CP 5 /ib) satisfies all experimental and observational bounds. Some of the remaining parame- 
ter space may be probed with forthcoming more sensitive direct searches and with the Gamma- Ray 
Large Area Space Telescope. 
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I. INTRODUCTION 



A wealth of observational evidence indicates the exis- 
tence of considerably more mass in galaxies and clusters 
of galaxies than we see in stars and gas. The source of 
the missing mass has been a problem since Zwicky's 1937 
measurement of the masses of extragalactic systems |l| . 
Given the evidence from galaxy clusters, galaxy dynamics 
and structure formation, big-bang nucleosynthesis, and 
the cosmic microwave background that baryons can only 
account for ~ 1/6 of this matter, most of it must be non- 
baryonic. Although neutrinos provide the cosmological 
density of dark matter if their masses sum to ~ 12 eV, 
such particles cannot (essentially from the Pauli princi- 
ple) have a sufficiently high phase-space density to ac- 
count for galactic dark-matter halos .2]; moreover, such 
masses are now inconsistent with neutrino-mass measure- 
ments |3J. Theorists have thus taken to considering for 
dark-matter candidates new physics beyond the standard 
model. To date, the most promising candidates — those 
that appear in fairly minimal extensions of the standard 
model and which coincidcntally have a cosmological den- 
sity near the critical density — are a weakly-interacting 
massive particle (WIMP), such as the neutralino, the su- 
persymmetricpartner of the photon, Z° boson, and/or 
Higgs boson U, or the axion pj. A considerable theo- 
retical literature on the properties and phenomenology 
of these particles has arisen, and there are considerable 
ongoing experimental efforts to detect these particles. 

In the absence of discovery of such particles, it may be 
well worth exploring other possibilities. Thus, an alterna- 
tive line of investigation takes a more model-independent 
approach and seeks to explore phenomenologically the 
possible properties of a dark-matter particle. Along these 
lines, for example, constraints to strongly-interacting 
dark matter were considered in Ref. Q; self- interacting 
dark matter has been considered Q , and some have stud- 
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ied whether dark matter might be charged Q or have a 
millicharge |9lll0|. 

Our investigation follows in spirit the latter possibil- 
ity. In particular, dark matter is so called because the 
coupling to photons is assumed to be nonexistent or very 
weak, or else we would have presumably seen such par- 
ticles either through absorption or emission of radiation 
or in laboratory experiments. In this paper, we ask the 
question, "How dark is 'dark'?" In other words, how 
weak must the coupling of the dark- matter particle to the 
photon be in order to be consistent with laboratory and 
astrophysical constraints? In the work on millicharged 
particles, a dark-matter coupling to photons was assumed 
to arise from a tiny charge. 

In this paper we consider the possibility that the dark 
matter possesses an electric or magnetic dipole moment. 
The punch line, illustrated in Fig. ^ is that a Dirac par- 
ticle with an electric or magnetic dipole moment of order 
~ 10~ 17 e cm with a mass between an MeV and a few 
GeV can provide the dark matter while satisfying all ex- 
perimental and observational constraints. 1 

In the following Section, we introduce the effective La- 
grangian for the dipolar dark matter (DDM) interac- 
tion with photons. We discuss the relic abundance in 
Section IIIII Section IIVI presents constraints on dark- 
matter dipole moments and masses that arise from di- 
rect searches at low-background detectors as well as con- 
straints from high-altitude experiments. Section [V] dis- 
cusses constraints due to precision tests of the standard 
model, while Section Ivfl discusses constraints due to the 
cosmic microwave background and the growth of large- 
scale structure. We provide some concluding remarks in 
Section IVIIII An Appendix provides details of the cal- 



1 Throughout, we will quote numbers for both the electric and 
magnetic dipole moments in units of e cm, where e is the electron 
charge. For reference, the Bohr magneton fig = eh/2m e = 
1.93 X 10 -11 e cm in these units. Also note that we work in 
rationalized Gaussian units so that the fine-structure constant 
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e /4ttTlc sa 1/137, and in particle-physics units (h = c = 1) 
a 4tt/137 and e pa 0.303. 
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FIG. 1: The constraints on the dipolar-dark-matter parameter space [m x , (D, M)\ that come from present-day searches and 
experiments. Viable candidates must lie in the shaded region, below the solid lines and outside the long-dashed lines. The 
short-dashed "relic abundance" curve shows where the dark matter would have a cosmological density Q x h 2 = 0.135, assuming 
standard freezeout, no particle-antiparticle asymmetry, and no interactions with standard-model particles apart from the dipole 
coupling to photons. 



culation of the drag force between the baryon and DDM 
fluids used in Section IVll 



II. THEORY OF DIPOLE MOMENTS 

A particle with a permanent electric and/or magnetic 
dipole moment must have a nonzero spin; we thus con- 
sider spin- 1/2 particles. Moreover, Majorana particles 
cannot have permanent dipole moments, so we consider 
Dirac fermions. Since the spin and the magnetic dipole 
are both axial vectors, a magnetic dipole moment can 
arise without violating any discrete symmetries. How- 
ever, the electric dipole moment is a vector and thus re- 
quires time-reversal and parity violation. 

The effective Lagrangian for coupling of a Dirac 
fcrmion \ with a magnetic dipole moment M. and an 
electric dipole moment V to the electromagnetic field F^ u 
is 
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At energies low compared to the dark-matter mass, the 
photon is blind to distinctions between M. and T> (un- 
less time- reversal- violating observables are considered). 



Hence, we can discuss limits to T> which equally apply to 
A4, except where noted. 

On dimensional grounds, we expect the electric dipole 
moment to be V < em" 1 ~ 2 x 10~ 14 (m p /m x ) e cm, 
where m p is the proton mass. Similar arguments also ap- 
ply to the magnetic dipole moment. This limit is shown 
as the perturbative bound in Fig. fj as violation of this 
bound would signal some non-trivial or non-perturbative 
field configuration. As we will see below, more rigorous 
but slightly weaker upper limits can be set with unitarity 
arguments. 

These upper limits already simplify our analysis. The 
phenomenology of charged dark-matter particles is de- 
termined largely by the ability of these particles to form 
atom-like bound states with electrons, protons, or each 
other. However, dipolar dark matter cannot form such 
bound states. A neutral particle with a magnetic moment 
will not form bound states with charged particles. Cu- 
riously enough, a neutral particle with an electric dipole 
moment (EDM) can indeed form a bound state with 
an electron, as first noted by no less than Fermi and 
Teller [111 > but only if the dipole moment is greater than 
0.639 e ao = 3.4 x 10~ 9 e cm (assuming m x 3> m e ) where 
ao is the Bohr radius. For smaller values of the dipole, the 
electron "sees" both poles of the dipole and finds no sta- 



ble orbit. This critical electric dipole moment scales in- 
versely with the dipole-electron reduced mass, so a bound 
state with a proton can occur if the dipole mass is ^> m p 
and T> > 1.8 x 10 _12 e cm. As we will see below, such 
values for the EDM cannot occur for a pointlike DDM. 
Likewise, the weakness of the dipole-dipole interaction 
prevents the formation of any stable dark-matter atoms. 
The first cosmological constraint is that from big-bang 
nucleosynthesis (BBN). BBN requires that the effective 
number of relativistic degrees of freedom at T ~ MeV 
does not exceed the equivalent of roughly 0.2 neutrino 
species [l2|. Since the particles we are considering 
are Dirac particles, they act like two effective neutrino 
species and thus cannot be relativistic and in equilib- 
rium at BBN. Generally, such considerations rule out 
m x < MeV, and so we restrict our attention in Fig. Q] 
to masses above an MeV. Strictly speaking, if the dipole 
moment is [D, A4] < 10~ 22 e cm, and if the particle has 
no other interactions with standard-model particles, then 
a particle of mass < MeV can decouple at a temperature 
> 10 GeV, and if so, it will evade the BBN limit. 



III. DARK MATTER ANNIHILATION AND 
RELIC ABUNDANCE 



DDM particles can exist in thermal equilibrium in the 
early Universe when the temperature T 3> m x , and their 
interactions will freeze out when T drops below m x re- 
sulting in some cosmological relic abundance. The mass 
density of relic DDM particles is fixed by the cross sec- 
tion er ann for annihilation to all lighter particles times the 
relative velocity v through (see, e.g., Eq. (5.47) in Ref. 
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where 



A = 0. 038 y/g^m p im x (cr ann v) 
6.3 x 10 9 (s*/!0) 1/2 



(m x /GeV) 



5.3 x 10 _26 cm 3 sec -1 
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assuming that annihilation takes place (as it does in our 
case) through the s wave. Here, g* is the effective num- 
ber of relativistic degrees of freedom at the temperature 
Tf ~ m x /A of freezeout. For the interaction of Eq. QJ, 
DDM-anti-DDM pairs can annihilate to either two pho- 
tons or to charged particle-antiparticle pairs through the 
diagrams shown in Figs. |21and|3J The cross sections for 
these two processes (to lowest order in v) are 

<r xx ^ 2l v = (X> 4 + M 4 )m 2 x /2tt 



a xx^ff v = N cB aCD 2 +M 2 ) 

= 2 x 10" 27 iV off (X> 2 7 + M\ 7 ) cm 3 sec" 1 , (4) 

where [T>n, M u ] = [V, X]/(10" 17 e cm), and m GcV = 
m x /GeV. Here, N c g — J2fQ 2 fN c f is the effective num- 
ber of fermion-antifermion pairs with mass to/ < m x , 
Qf is the charge of fermion /, and N c f is the num- 
ber of color degrees of freedom for fermion /. (N c f = 
1 for electrons.) In the standard model, annihilation 
can also occur to W + W~ pairs above threshold. For 
(2?i7, Adn)(m x /m p ) < 5000, fermions are the dominant 
final state. The present-day mass density of DDM par- 
ticles thus depends primarily on the dipole moment. If 
such particles are to account for the dark matter, then 



n x h 2 = 0.135 0, and (V 2 + M 2 ) 1 ' 2 



1.0 xlO" 17 , 



for TOv 



1 GeV. The full mass dependence of this result 



is shown in Fig. ^ 

The cross sections in Eq. Q are s-wave cross sections. 
According to partial-wave unitarity, the total s-wave an- 
nihilation cross section must be a < 47r/m 2 [l5[ , so that 
(T>, Ai)m x < 3, fixed by the cross section for annihilation 
to two photons. This limit is shown in Fig. Q] as is the 
more stringent, but less rigorous, limit (T>, M) % e/m x . 

Of course, the present-day mass density of DDM par- 
ticles could differ from the estimates obtained above. To 
obtain these results, we assumed (1) that the dipole in- 
teraction with photons is the only interaction of these 
particles; and (2) that there is no particle-antiparticle 
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FIG. 2: Feynman diagrams for annihilation of a DDM-anti- 
DDM pair to two photons. 




FIG. 3: Feynman diagram for DDM-anti-DDM annihilation 
to fermion-antifermion pairs. 



asymmetry. It is reasonable to assume that in any realis- 
tic model, a dark-matter dipole interaction will arise from 
loop diagrams involving other standard-model and new 
particles. If so, then there may be other contributions 
to the annihilation cross sections. In this case, the relic 
abundance will be smaller than we have estimated above. 
We thus conclude that if there is no particle-antiparticle 
asymmetry, our estimates should be treated as an upper 
bound to the relic abundance, and the Q, x h 2 curve in Fig. 
n should thus be considered an upper limit to the desired 
values of T> and Ai . On the other hand, the relic abun- 
dance could also be increased if exotic processes increase 
the expansion rate during freezeout jlqj . 

If there is an asymmetry between \ an d X, then the 
relic abundance is increased relative to our estimate. In 
this case, however, the present-day Universe should con- 
tain predominantly either particles or antiparticles. Al- 
though there is no a priori reason to expect there to be 
a particle-antiparticle asymmetry, the observed baryon- 
antibaryon asymmetry might lead us to expect an analo- 
gous dark-matter asymmetry, should the dark matter be 
composed of Dirac particles. It is possible such asymme- 
tries have a common origin. 

Finally, we have assumed above that the particles 
freeze out when they are nonrelativistic. However, as 
the dipole moment is lowered for a given mass, freezeout 
occurs earlier. If the dipole moment is reduced beyond 
a certain value, and if there are no other couplings to 
standard-model particles, then freezeout will be relativis- 
tic. These particles will then be roughly as abundant as 
photons, and they will overdose the Universe by huge 
margins unless their masses are < few eV; even in this 
case they will violate constraints to hot dark matter from 
the CMB and large-scale structure, and they will also be 
unable, from the Tremaine-Gunn argument, to make up 
the dark matter in Galactic halos. The transition from 
nonrelativistic to relativistic freezeout occurs (again, as- 
suming no non-dipole interactions with standard-model 
particles) for m x T>\ 7 < 10~ 10 GeV for m x > MeV, and 
for m x < MeV, at m x V% Z < 200 MeV. 



IV. DIRECT DETECTION 

The diagram for scattering of a DDM particle with 
a particle of charge Ze occurs through the exchange of 
a photon, as shown in Fig. 0] (not unlike the electron- 



neutron interaction 17]). In the nonrelativistic limit, the 
differential cross section for this process is given by, 
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FIG. 4: Feynman diagrams for elastic scattering of an electron 
from a DDM particle. 
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where v is the relative velocity. Through this interaction 
the DDM might produce a signal in a direct-detection 
experiment. Although the total cross section is formally 
infinite, the divergence comes from the small-momentum- 
transfer scatterings that will be screened by atomic 
electrons. Roughly speaking, then, the DDM-nucleus 
elastic-scattering cross section will be a ~ (Ze) 2 (T> 2 + 
M 2 )/2nv 2 ~ 6.4 x 10~ 32 Z 2 {V\ 7 + M\ 7 ) cm 2 , using 
v ~ 10~ 3 c. Current null searches in germanium de- 
tectors [(A, Z) = (76, 32)] correspond for masses m x ~ 
10 GeV to a rough upper limit to the cross section 
~ 10 -42 cm 2 llq . thus ruling out any dipole moment 
(V\ 7 + M 2 17 ) l/2 > 10~ 7 . This is shown in Fig. Q] as 
the horizontal dashed line at T> = 10 _24 e cm. Note that 
the cross-section limit depends (and increases) with mass 
at higher masses; the curve appears as a horizontal line 
simply because of the break in scale on the y axis. 

This seems like a very stringent limit, especially consid- 
ering the value, T>\ 7 ~ 1, favored for the correct cosmo- 
logical density. However, if the dipole moment becomes 
large enough, the particles will be slowed in the rock 
above the detector and thus evade detection in these un- 
derground experiments. In order to determine the mag- 
nitude of the dipole moment for this to occur, we need to 
calculate the stopping power dE/dx for the particle as it 
passes through the atmosphere and then the rock. Since 
elastic scattering takes place through exchange of a pho- 
ton, it leads to a long-range interaction and, as we have 
seen above, a formally divergent elastic-scattering cross 
section. The calculation of the stopping power thus par- 
allels that for ionization loss due to Coulomb collisions, 
with two important differences. First of all, since the 
long-range force is oc r -3 , as opposed to r -2 for Ruther- 
ford scattering, stopping occurs via scattering from nu- 
clei, rather than electrons. Second, since this interaction 
falls more rapidly with radius than the Coulomb interac- 
tion, the stopping power is due primarily to hard scatters 
at small impact parameter, rather than soft scatters at a 
wide range of impact parameters. 

Our result for the stopping power due to scattering 
from nuclei of charge Z is 



dE 
dx 



-n N 



(Ze) 2 (V 2 +M 2 Wc 2 

Tda = -n N , (6) 
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where the kinetic-energy transfer in a single collision 
is T — p (1 — cos 9) /rriN, x is the depth, and /i = 
//[m x ,mjv] = m x mN(rn x + m_/v) _1 is the reduced mass. 
For very weak WIMP interactions with nuclei, the most 
restrictive limits on the WIMP-proton cross section (the 
smallest upper bounds) are obtained from null searches 
from experiments that are deepest underground (so as 
to reduce the background). However, the most restric- 
tive constraints on the cross section at the upper end of 



the excluded range of cross sections will come from the 
shallowest underground experiment with a null result. 
From Eq. ©, we find that the dark-matter particles will 
only penetrate to a depth x = (E{ — Ef)/\dE/dx\ where 
Ei = 2 m x y2 i s the initial dark-matter kinetic energy and 
Ef is the final energy. For a stopped particle, Ef = 0. 
However, the particle only needs to lose enough energy 
for Ef to drop below the detection threshold for a partic- 
ular experiment. Equating the maximum kinetic energy 
transferable in a collision (8 — tt) to the threshold de- 
tectable nuclear-recoil energy (-Eth), we find the velocity 
must be slowed to vi = mdEt\j2^[m X) md] 2 , where rrid 
is the mass of the nuclei in the detector, and fi[m x ,m,d] is 
the DDM-nucleus reduced mass. Hence, the final dark- 
matter energy must be Ef = m x md-Eth/4/i[m x , m^] 2 . 

So far we have assumed that the particle loses en- 
ergy but is not significantly deflected in each scatter; 
this will be a good approximation if m x 3> tojv- How- 
ever, when m x < msi — 26 GeV, the dark-matter par- 
ticle may be backscattered upon encountering a terres- 
trial nucleus, rather than simply being slowed without 
deflection. In this case, the particle will diffuse, under- 
going N ~ mjv/2m x scatters before coming to rest. If 
so, the penetration depth will be reduced by an addi- 
tional factor of ~ TV -1 / 2 . We thus replace dE/dx —* 
dE/dx[l + (mN /2m x ) 1 ' 2 ] in our expression for the pen- 
etration depth. 

Then, for a given shielding thickness L, in meters wa- 
ter equivalent (mwe), we invert the expression for the 
stopping distance to obtain the following bound on the 
dipole strength: 



V 2 +M 2 > 
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where the index i sums over the composition of the shield- 
ing material, and fi is the fractional composition by 
weight. We use a realistic model of the composition of 
the Earth (chemical composition by weight: O [46.6%], 
Si [27.7%], Al [8.1%], Fe [5%], Ca [3.6%], Na [2.8%], K 
[2.6%], Mg [2.1%]) and atmosphere (10 mwe consisting of 
a 4:1 ratio of nitrogen to oxygen), although the resulting 
bounds do not change substantially if we ignore the at- 
mosphere and approximate the Earth's crust as entirely 
composed of Si. We take the initial DDM velocity to be 
300 km sec -1 . 

The shallowest underground experiment with a strong 
null result is the Stanford Underground Facility (SUF) 
run of the Cryogenic Dark Matter Search (CDMS) [l9| . 
which was situated at a depth of 16 mwe. With a de- 
tector energy threshold of E t h = 5 keV, it is sensitive to 
DDM masses down to m x ~ 10 GeV. Near this threshold 
we find that DDM particles are stopped by the shielding 
for 2?i7 > 20. This bound grows more prohibitive with 
increasing mass, as illustrated in Fig. Q] The Cryogenic 
Rare Event Search with Superconducting Thermometers 
(CRESST) 20], though at a depth of 3800 mwe, extends 
to slightly lower masses, having a detector threshold en- 



ergy Eth = 0.6 keV. Near m x ~ 1 GeV the minimum 
dipole strength is T>n > 2. However, there are no limits 
from underground experiments for DDM masses below 
1 GeV. 

Two airborne experiments — unobscured by the atmo- 
sphere or rock — which have closed the windows on some 
forms of strongly- interacting dark matter [g, l2lj , also re- 
strict dark-matter dipole moments. To determine the 
predicted signal at a detector, we recast Eq. J5j as 
the cross section per energy transfer, whereby da/dT = 
Z 2 e 2 (V 2 +M 2 )/4ttv 2 T. The event rate (per time, en- 
ergy, and unit mass of detector) is 



R = 7X^(0.3 GeV cm" 



v da 
m x dT 
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keV 
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where iVjy is the number of nuclei per gram of material. 

The silicon semiconductor detector flown on a balloon 
in the upper atmosphere by Rich et al. | 22| observed an 
event rate of ~ 0.5 counts sec -1 keV -1 g _1 nuclear recoils 
in the lowest energy bin at 2 keV. For dark-matter masses 
above the threshold ~ 7 GeV, we thus require (V\ 7 + 
M 2 17 )(m p /m x ) < 0.2. 

The X-ray Quantum Calorimeter (XQC) detector 
flown on a rocket by McCammon et al. J23J, was de- 
signed to probe the soft x-ray background. However, 
it serendipitously provides a tight constraint on dipo- 
lar dark matter. To predict the expected number of 
events, we start by computing the number of DDM parti- 
cles that could impinge on the detector: N x = n x vAt = 
3 x 10 7 m p /m x , where n x is the galactic number density 
of dark-matter particles, v is their velocity, the cross- 
sectional area of the XQC detector is A = 0.33 cm 2 , 
and the rocket flight was about t ~ 120 seconds. The 
chief property of the XQC detector is the 14 //m thick 
Si substrate above the thermistor, providing a target of 
-/Vjv ~ 6.5 x 10 19 nuclei/cm 2 . Thus, the event/energy 
count NNN x da/dT integrated over the 25-100 eV en- 
ergy bin gives a predicted ~ 0.38 (D 2 7 + M\ 7 )(m p /m x ) 
events compared to the ~ 10 observed events. Since the 
detector has a 25 eV threshold, energy transfer by dark- 
matter particles as light as ~ 1 GeV can be detected. 
Altogether, the balloon and rocket experiments exclude 
a wide range of dipole strengths and masses, as illustrated 
in Fig. U| 



V. CONSTRAINTS FROM PRECISION 

MEASUREMENTS 



We now consider the limits placed on DDM due to 
precision tests of the standard model. Our use of per- 
turbation theory is valid provided the energy scale of the 
interaction £ satisfies (T>, Ai)£ -C 1. In addition, we 
require that the DDM mass satisfies (T>, Ai)m x < 1, 
equivalent to the unitary bound |15|. which ensures the 
self-consistency of the local operator in Eq. QJ . Indeed, 



if A is the energy scale at which a dipole is generated then 
one generally expects (T>, A4) A ~ 1. In £ 7X we assume 
that all interacting particles with masses greater than A 
have been integrated out. Consequently, one must have 
at least m x < A for the dark matter to be dynamical, 
which also yields (D, M)m x < 1. 



A. Muon Anomalous Magnetic Moment 

The interaction in Eq. Q contributes to the pho- 
ton propagator via the diagram shown in Fig. [5] The 
photon-DDM interaction vertices are either both electric 
or magnetic dipolar; the mixed diagram where one ver- 
tex is magnetic and the other is electric is proportional to 
^^up\F^ v ' F pX — for photons with equal momenta. The 
sum of the diagrams produces the following contribution 
to the photon vacuum-polarization tensor: 



W{q 2 ) 



(3 = 



n(<? 2 ) (^v 

V 2 +M 2 ( 
8tt 2 



1, mi' 
1 - r In — £ 



fiq 2 (g V - qVq") 



(9) 



where the photon momentum is taken to be small, q 2 <C 
m 2 (resulting in j3q 2 <1), and /i is the renormalization 
scale, which should be smaller than A. We take \x < 1 
TeV for our estimates. With this self-energy correction, 
the photon propagator for (3q 2 <C 1 can be written as 



■'""" '? 



1 



1//3 



(10) 



The second term above generates a correction to the 
muon gyromagnetic ratio Sa^ = —am 2 /3/3ir. Interest- 
ingly, this contribution is not explicitly suppressed by 
the DDM mass. In view of recent measurements J24J 
and comparison with the SM predictions, we require that 
Sa^ does not exceed 10 -9 , whereby (V 2 + M 2 ) 1 ^ 2 < 
6 x 10~ 15 e cm. The order of magnitude of this result can 
be obtained on dimensional grounds, if we consider that 
the DDM dipole moment contributes to a M via at least 
a two- loop graph (see Fig. 0) , with two electromagnetic 
couplings and two dipole couplings. Including a factor 
1/167T 2 per loop, one obtains the estimate, 



Sa L 



(16tt 2 ) ; 



(V 2 +M 2 )£ 2 , 



(11) 



where £ is the characteristic energy scale for the process. 



In the case of the muon, £ 



1 1 > 



which reproduces the 



rigorous result to within an order of magnitude. 



B. Electric Dipole Moments 




FIG. 5: One-loop correction to the photon self-energy induced 
by dipole moments A4,T> of the dark-matter particle. 

X 




FIG. 6: Lowest-order correction to the muon anomalous mag- 
netic moment induced by dipole moments of the dark-matter 
particle. 



closed loop. 2 Contributions with an odd number of pho- 
tons attached sum to zero. On the other hand, one must 
have an odd number of time-reversal-odd (T-odd) EDM 
vertices in the DDM loop to generate a T-odd opera- 
tor. These considerations show that the lowest-order 
non- vanishing diagram must have four photons attached 
to the DDM loop; diagrams with two photons attached, 
similar to the one in Fig. vanish (see above). Out 
of the four photons attached to the DDM loop, either 
one or three can have EDM coupling to DDM. Note that 
in this scenario both electric and magnetic DDM mo- 
ments are necessary to generate a dipole moment for 
a SM fermion. With these considerations in mind, the 
lowest-order three-loop diagram that induces an EDM 
for a charged fermion is shown in Fig. One obtains 
the following estimate for the induced electric dipole mo- 
ment: 



X>/~ [DM{V 2 +M 2 )] 



e 3 m 3 



/ 



(16tt 2 ) 



2\3 



In 2 



mj 



(12) 



where a possible double-logarithmic enhancement is in- 
cluded. For the electron, the present limit is T> e < 
4 x 1CT 27 e cm, which implies (V, M) < 3 x 10~ 13 e cm 



for m x ~ 100 GeV and V 



M.. For smaller m x the 



limit becomes weaker. 

There are constraints on the EDMs of other systems, 
such as the neutron and the mercury atom. It is non- 
trivial to translate such constraints into limits on the un- 
derlying interaction. In case of the neutron, one may at- 
tempt to treat the neutron as a point particle for virtual- 
photon energies below 1 GeV. For higher loop momenta, 



Furry's theorem tells us that in evaluating radiative 
corrections to a process one should only keep the dia- 
grams with an even number of photons attached to a 



2 Since the theorem is valid for interactions that preserve charge 
conjugation invariance we can apply it to electric and magnetic 
dipole moment interactions. 




FIG. 7: Three-loop contributions to the EDM of a charged 
fermion /. Either one or three of the DDM-photon interac- 
tion vertices must be of EDM type. The dots indicate all 
other diagrams which can be obtained from the one shown by 
permutation of the interaction vertices. 




FIG. 8: Lowest-order correction to Z°-pole observables in- 
duced by dipole moments of the dark-matter particle. 




FIG. 9: The photon-DDM scattering diagram, the analog of 
the Compton-scattering diagram for electric or magnetic mo- 
ments. 



pression for Ar, whereby Ar 
the standard model A: 



New 



,SM _ 



n(m|) - 11(0). In 
0.0355 ±0.0019 ±0.0002. 
On the other hand, one can use experimentally mea- 
sured values for a, mw,z, and Gf in Eq. (|14|) to infer 
Ar Ex P = 0.0326 ±0.0023, which gives Ar New < 0.003 at 
the 95% confidence level. Therefore, we obtain the limit 
(V 2 + M 2 ) 1/2 < 3 x 10" 16 e cm. A full calculation of 
the vacuum polarization yields the constraint shown in 
Fig. n This turns out to be the strongest constraint due 
to precision tests of the standard model. 



photons begin seeing the quarks and the contribution to 
the EDM becomes suppressed by the quark masses. In 
this case one may use the above equation with rrij — > m n 
and no logarithmic enhancement, in order to estimate the 
neutron EDM: 



V n = [VM(V 2 +M 2 )] 



e \Kt, 



| 3 m 3 



(16tt 



2\:s 



- < 6 x 10~ 26 e cm, 



(13) 



which results in the limit (T>, JVl) < 4 X 10 -15 e cm (as- 
suming T> ~ M.. In the above equation, n n = —1.91 
is the magnetic moment of the neutron. It appears be- 
cause the neutron is neutral, and couples to the photon 
in Fig. [7] via a magnetic-moment interaction. The limit 
for the EDM of the mercury atom is much stronger than 
the neutron, T>n g ^$ 10~ 28 e cm. Unfortunately, the mer- 
cury atom is a complicated system for which the EDM is 
influenced by many sources. Therefore, we leave the Hg 
limit for future study. 



C. W Boson Mass 

The DDM can contribute to the running of the fine- 
structure constant for momenta ranging up to the Z° 
mass. Such running will affect the relationship between 
the Fermi constant Gf, the mass of the W^ boson, and 
the fine-structure constant at zero momentum: 



1 



l w 



V2G F (1 " <K)(1 - Ar) : 



(14) 



where Ar is a correction calculable in a given theory. 
The interaction in Eq. JJ| modifies the standard ex- 



D. Z-Pole Observables 

The DDM will contribute to various Z°-pole observ- 
ables through two-loop diagrams similar to the one shown 
in Fig. El at the level a(D 2 + M 2 )m 2 z /64Tr 3 . Requiring 
that these contributions do not exceed the ~ 0.1% pre- 
cision to which Z°-pole observables are typically known 
[1| results in the constraint (V 2 +M 2 ) 1/2 < 10~ 14 e cm. 
Note that in order for perturbation theory to apply for 
energies £ ~ mz, one must have (T>, A4)mz < 1, which 
means (V, M) < 7 x 10 -16 e cm. Interestingly, con- 
sistency with a perturbative treatment at the Z° pole 
imposes much stronger constraints on the DDM than the 
Z°-pole observables themselves. 



E. Direct Production 

If kinematically allowed, DDM can be directly pro- 
duced in various scattering and decay experiments. In 
this case one may use the "missing energy" signature to 
constrain the DDM couplings. Here, we consider missing- 
energy constraints from both low-energy (B + and K + 
meson decays) as well as collider (LEP, CDF) experi- 
ments. 



B+ and K+ 



Searching for light (rn x < 1 GeV) dark matter us- 
ing missing-energy signatures in rare B + meson de- 
cays was originally suggested in Ref. J2|j. There, data 
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FIG. 10: 
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collectively represent quark flavor-changing interactions. 




v,X 



v,X 



from BABAR J2g and CLEO |2g were used to set a 
limit Br(B+ — > iiT + + invisible) < 10 -4 [derived from 
Br{B + — ► iif+i/j/)]. This limit can be used to constrain 
the dipole moments of dark matter. The diagram for 
_B + — » iiT+xx decay is shown in Fig. llUf aL The rate for 
this decay can be related to the photon-exchange contri- 
bution to B+ -> ir+/+/"| 7 shown in Fig. E^b). Since the 
graphs have identical topologies, the difference in rates 
will come from the difference in effective couplings and 
the final-state phase-space integrals. One can estimate, 



Br (B+ -> K+xx) 
Br{B+->K+l+l-)\~ 



Dffli 



PS (K+xx) 



ps(K+i+i-y 

(15) 
where PS(- ■ ■ ) stand for the corresponding final-state 
phase-space integrals, and m B + = 5.279 GeV is the 
B+ mass. Belle [29J and BABAR H3 find Br(B+ -> 
K+l+l~) < 10~ 6 . Since the ratio of the phase-space in- 
tegrals is of order unity, and since in the absence of acci- 
dental cancellations Br(B+ -> K+l+l-% < Br(B+ -v 
K+l+l~), one obtains the constraint 



2 x i.r 6 f'^^i ) < 10- 4 



(16) 



which leads to I? < 3.8 x 10~ 14 e cm. This constraint is 
relevant for m x < (rriB+ — m^+)/2 = 2.38 GeV. 

Rare K+ decays can be treated in a similar manner. 
The relevant branching ratios are Br(K+ — > 7r + e + e~) = 
2.88t^a x iO" 7 and Br(if+ -> tt+Pi/) = 0.157^82 x 
10 -9 J25| . The resulting constraint on the dipole moment 
is V < 1.5 x 10 -15 e cm. This constraint applies for m x < 
[mx+ — "t-tt+)/2 = 0.18 GeV. We see that constraints 
from B+ and K+ decays are not competitive with other 
constraints shown in Fig. E 



2. Collider experiments 

A typical example of a process where DDM can be 
directly produced in a collider experiment is shown in 
Fig. ^2 Here, two fermions / scatter to produce a final 
state containing some set of visible particles X (photon, 
multiple jets, etc.) along with particles that are not de- 
tected. In the SM, the latter are neutrinos. Limits on 



FIG. 11: A typical process that would produce a missing- 
energy signature in a collider experiment. Here, X stands for 
the visible portion of the final state. Neutrinos or DDM may 
carry a large fraction of the energy but are not detected. 



the rate for such processes have been set by, e.g., the L3 
and CDF collaborations 31]. At LEP, X consisted of a 
single photon whereas at CDF it consisted of one or more 
hadronic jets. 

In order to translate constraints from collider exper- 
iments into limits on DDM couplings one needs an an- 
alytical expression for the rate for // — ► Xxx- Naive 
application of the effective Lagrangian in Eq. (JJJ would 
result in upper limits from these missing-energy searches 
of roughly 10~ 17 e cm. However, this constraint does not 
actually exclude larger values of the dipole moments. In- 
deed, as discussed above, perturbation theory will break 
down when the energy scale for the process £ satisfies 
T>£ > 1. This means that missing-energy searches from 
L3 {£ w 200 GeV) and CDF {£ = 1.8 TeV) cannot be 
used to probe effective dipole moments V > 10~ 16 e cm 
and T> > 10~ 17 e cm, respectively, unless the underlying 
physics that gives rise to the dipole moment is specified. 



F. Other Laboratory Constraints 



Important constraints can be obtained for millicharged 
particles from the Lamb shift j^, yj] and from a tar- 
geted experiment at SLAC J3T|. We have checked, how- 
ever, that due to the different energy dependence of 
the photon-dipolc vertex, as opposed to the photon- 
millicharge vertex, the DDM-induced correction to the 
Lamb shift is small for dipole moments not eliminated by 
other precision measurements, such as the running of the 
fine-structure constant. Likewise, although the SLAC ex- 
periment is in principle sensitive to neutral particles with 
a dipole, the production and energy deposition of dipole 
particles is sufficiently small, for dipole moments consis- 
tent with accelerator experiments, to evade detection in 
the SLAC experiment. 



VI. CONSTRAINTS FROM LARGE-SCALE 
STRUCTURE AND THE CMB 

We now consider the effects of the interaction L lx on 
the evolution of cosmological perturbations and their re- 
sulting imprints on the matter power spectrum and the 
CMB. A dipole moment can induce a coupling of the 
dark matter to the baryon-photon fluid by scattering 
from photons through the diagrams shown in Fig. |5J or 
by scattering from protons, helium nuclei, and/or elec- 
trons through the diagram shown in Fig. 21 What we 
will show below is that the dark matter is coupled to 
the baryon-photon fluid at early times, and decouples at 
later times. When the dark matter is coupled to the 
photon-baryon fluid, the pressure of the plasma resists 
the growth of gravitational potential wells. Thus, the 
short-wavelength modes of the density field that enter 
the horizon at early times will have their growth sup- 
pressed relative to the standard calculation resulting in 
a suppression of small-scale power. The evolution of the 
longer- wavelength modes that enter the horizon after the 
dark matter has decoupled remain unaffected. Before 
presenting the results of our detailed analysis, we begin 
with some simple estimates. 

We first show that DDM-photon scattering is negligi- 
ble compared with DDM-baryon scattering in providing 
the drag force between the DDM fluid and the baryon- 
photon fluid. To do so, we first estimate the drag force 
per unit mass (i.e., the deceleration) on a DDM particle 
that moves with a velocity V with respect to the rest 
frame of a blackbody at temperature T. The diagrams 
in Fig. will lead to a photon-DDM scattering cross 
section a xl ~ T> A E 2 Considering that the momentum 
transfer to the DDM particle in each scatter is ~ E 1 and 
that the difference of the fluxes of photons moving in 
the same versus opposite direction to the DDM particle 
is ~ T 3 (V/c), we conclude that the deceleration due to 
photon scattering is a 7 ~ T> 4 T 6 (V/c)/m x . 

We next estimate a&, the drag force per unit mass due 
to DDM-proton scattering. We first note that the pe- 
culiar velocity of the baryon-photon fluid (obtained from 
the continuity equation) in the early Universe will be V ~ 
(H/k)cS, where H is the Hubble parameter, k is the phys- 
ical wavenumber of the mode in question, and S ~ 10 
is the amplitude of the fractional density perturbation. 
Since {H/k) < 1 for modes inside the horizon, we must 
have V < 10 -5 c. On the other hand, the proton thermal 
velocity dispersion is v p ~ (T/m p ) l l 2 c > 10~ 45 c before 
recombination. Thus, for the early times of interest to 
us here, the relative velocity between the DDM and the 
baryon-photon fluid is small compared with the thermal 
proton velocities. Thus, the appropriate relative velocity 
to use in Eq. JSJ in estimating the proton-DDM cross 
section is v p , resulting in a DDM-proton cross section 
<7 X 6 ~ e 2 T) 2 /v 2 The momentum transfer per scatter is 
~ /j,v p , where /z is the proton-DDM reduced mass, and 
the difference of the fluxes of protons moving with as 
opposed to against the DDM fluid is n p V, where n p is 



the proton density. The drag force per unit mass on 
the DDM fluid due to scattering with protons is thus 
a& ~ e 2 V 2 (fi/m x )(V/v p )n p . We also conclude from the 
appearance of /i in this result that drag due to scattering 
from electrons is negligible compared with baryon drag. 
Since n p ex T 3 and v p ex T 1 ' 2 , we find a& ex T 25 as 
opposed to a 7 ex T 6 . Thus, at early times, photon drag 
dominates while baryon drag dominates at later times. 
The transition occurs at a temperature T ~ GeV for val- 
ues of m x and T> of interest to us, and such high temper- 
atures correspond to (comoving) horizon scales consid- 
erably smaller than the distance scales (> Mpc) probed 
by large-scale structure. We can thus neglect photon 
drag. From a& ex T 25 we infer a deceleration time for 
the DDM fluid t dcc = V/a b ex T~ 2b . Since this de- 
creases more rapidly than the Hubble time tu ~ mp\T~ 2 
(where mp\ ~ 10 19 GeV is the Planck mass), we conclude 
that DDM particles are tightly coupled to the plasma 
at early times and then are decoupled at later times. 
With these rough estimates, the transition temperature 
is T ~ 10 keVP 1 " 5 4 (l+m x /m p ) 2 suggesting that power on 
scales smaller than A ~ 10 -2 T>f 5 (l+m x /m p )~ 2 Mpc will 
be suppressed. The T 05 dependence of the ratio of the 
deceleration and expansion times suggests furthermore 
that the small-scale suppression will change gradually, 
rather than exponentially, with wavenumber k. Know- 
ing that the linear-theory power spectrum is measured 
and roughly consistent with scale invariance down to dis- 
tances ~Mpc leads us to conclude that dipole moments 
V > 5 x 10~ 15 (1 + m x /m p ) 1 / 2 e cm will be ruled out. 
Strictly speaking, when m x < m p , the detailed calcula- 
tion must take into account the velocity dispersion of the 
DDM particles; our detailed calculation below includes 
these effects. As seen below, the detailed analysis leads 
to a slightly stronger constraint. 



A. Exact Equations 

The standard calculation of perturbations in an ex- 
panding universe requires solution of the combined Ein- 
stein and Boltzmann equations for the distribution func- 
tions of the dark matter, baryons, photons, and neutrinos 
including all relevant standard-model interactions (see, 
e.g., Refs. HH, |35| and references therein). Since the 
perturbations are initially very small, linear perturba- 
tion theory is an excellent approximation; this allows us 
to solve the perturbation equations in Fourier space at 
each wavenumber k independently of all other wavenum- 
bers (modes are uncoupled). The scattering of photons 
and baryons by DDM through the interaction £ 7X influ- 
ences the growth of cosmological perturbations by intro- 
ducing additional collision terms to the Boltzmann equa- 
tions, which ultimately result in a drag force between the 
DDM and the colliding species in the equations describ- 
ing the cosmological fluid (see, e.g., Refs. Ha,|33, which 
consider similar effects) . Below we present the exact per- 
turbation equations including the effects of dark matter 



10 



with electric or magnetic dipole moments. Since solutions 
to these equations are numerically intensive when pho- 
tons and baryons are tightly coupled through Compton 
scattering, we also discuss the equations appropriate for 
solving for the DDM, photon, and baryon perturbations 



during the epoch of tight coupling. 

In the synchronous gauge the equations describing the 
evolution of baryons, photons, and dark matter with an 
electric or magnetic dipole moment are 



J 



2- 



— " - — h . dh = — Oh — — h . o v — — p v — — h . 

3 7 3 6 2 ' x x 2 ' 

k 2 I -<5 7 — £ 7 I + an e <jT(&b — 0y) + an x( a )xi^x ~ ^y)> 



9 b + c? sb k 2 5 b + -—^an e a T {0j - b ) + an x (av) xb (0 x - 6 b ), 



3/96 

Pb 

1 

Px 



\ -r c 2 sx k 2 5 x + ^-an x (av) xb (6 b - 6 X ) + ^an x (a) X7V , 



3Px 



(17) 



r 



While the evolution equations for the density contrast 
Sj = Spj/pj for each species j € {7,6, x} are as in 
the standard case [34], as discussed above, the evolu- 
tion equations for the fluid-velocity perturbations have 
additional drag-force terms due to the photon-DDM in- 
teraction. Note that in these equations and what follows 
the variable 8j = ikVj is the divergence of the fluid ve- 
locity in Fourier space, E-,- is the shear, c s j is the intrinsic 
sound speed, and rij and pj are the background number 
and energy densities of a particular species j, respec- 
tively. The variable h is the trace of the scalar metric 
perturbation in Fourier space (not to be confused with 
the Hubble parameter), a is the cosmological scale factor, 
and an overdot represents a derivative with respect to the 
conformal time r. Furthermore, <jt is the Thomson cross 
section, while 



1x1 



= ~^T 2 

21 7 ' 



(18) 



is the appropriately thermally-averaged DDM-photon 
cross section, which can be obtained from the differential 
cross section [38, 391, 



da 



_ (V 4 + M 4 )E 2 

dfT ~ 8^ 



(3 - cos 2 



(19) 



for photon-DDM scattering. As argued above, the 
photon-DDM drag term is small, and we consider it no 
further in Eq. itTTT) . 
The quantity 

; x 4(1 + £Y) m x o . , , J? . 
(crv) xb = ; g ' = * — e 2 (V 2 + M 2 ), 

,X 3ttV<^) 2 + K) 2 ™x+™p 

(20) 

is the appropriate thermally-averaged cross section times 

relative velocity for the baryon-DDM coupling, and 




m x + 4m p y {v p } 2 + 4(u x ) 2 



- 1 



(21) 



is the relative efficiency for coupling to helium nuclei 
compared to protons. Appendix A provides a derivation 
of this collision coefficient. In these expressions, Y — 
PHel Pb — 0.24 is the cosmological helium mass fraction 
(approximating mHe — 4m p ), (v p ) = \/8T b /irm p is the 
average thermal speed of the protons, (v x ) = \/8T x /iTm x 
is the average thermal speed of the DDM, and T 7 , T b , 
and T x are the photon, baryon, and DDM temperatures 



respectively, 
cording to 



The dark-matter temperature evolves ac- 



T x = 



a 2ap b {av)^ b 



-2-T, 



(T b - T x 



3m- 



^1(T 7 -T X ), 



(22) 



where (o"f)3b is the same as the expression given in Eq. 
(|20|l with the replacement of £ by £ T which is given by the 
expression in Eq. I|21|) with the factor (m x +m p )/(m x + 
4m p ) replaced by [(m x + m p )/(m x + 4m p )] 2 . The final 
term, describing the dark matter heating by photons, is 
important at very early times. For the dipole strength 
and mass range considered, the influence is manifest only 
on very small length scales, below the range of interest. 

At early times, the DDM temperature T x ~ T b , but 
at later times, when the DDM decouples, T x drops rel- 
ative to T b . The DDM-proton cross section is ex v~ 2 , 
which leads to (<Jv) xb ex (w p ) _1 . As a result, we can- 
not directly apply the results of Ref. [3||, wherein a 
velocity-independent dark matter-baryon interaction was 
assumed. However, we have verified that we recover their 
results if we take a velocity-independent cross section as 
the source of dark-matter-baryon drag. 
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B. Tightly Coupled Equations 



At early times when r" 1 = an e cr*r 3> a/a the rapid 
scattering of baryons and photons forces these species to 
have nearly equal fluid velocities, and consequently the 
solution of the equations shown in Eq. (|17(l is numeri- 
cally intensive. Following standard procedures J3J, |4(| 



we derive a set of equations to leading order in the (con- 
formal) Compton scattering time r c that are appropriate 
for evolving the fluid variables through this epoch of tight 
coupling. We first write down an equation for the time 
derivative of 6b — 6~ t which is usually termed the baryon- 
photon 'slip' to leading order in t c , 
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2R ib a,Q 



1 + R 



•yb ' 



1 + R 



*yb 



•V 

a 



k 2 f c 2 sb S b - -5^ 



r 



2a 1 



2or v l x 



,(23) 



where we have introduced the (conformal) DDM-baryon 



scattering time t x 



an x (<jv) x bi and i? 7 j, = 4p-y/(3pb) 



It is useful to separate this equation as a sum of the terms 
not containing t x (this is just the time derivative of the 
standard slip, which we denote SVy), and the new terms 
introduced by the DDM coupling, 



is the parameter that controls how strongly the new in- 
teraction affects the evolution of the slip. In terms of 
these definitions the baryon-velocity evolution equation 



= S, 



by 



•/? 



?x -Ob) - 2~(^x ~9b) 



, (24) 



where 
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(26) 



The photon-evolution equation is then given by the exact 
expression 



1 

Rfb 



h + -Ob-c 2 bk 2 S b -—(0 x -6b) 
a t v 



(27) 



k I -5-y — S 7 



We use these equations to follow the initial evolution of 
the baryon and photon fluid variables and switch to the 
exact equations of Eq. I|17(l at later times. For the evolu- 
tion of the DDM fluid variables we always use the exact 
form of Eq. (fTTjl . 



C. Effects on the Matter and CMB Power Spectra 

In Fig. 1121 we show the linear matter power spectrum 
and in Fig. 1131 we show the angular power spectrum of 
the CMB for several values of the dipole moment and 



for DDM mass m x = 1 GeV. Physically, the effects of 
DDM on the matter power spectrum and CMB can be 
simply understood. Prior to matter-radiation equality 
the photons have a much larger density than the baryons 
or the DDM and so to a first approximation completely 
drive the behavior of the baryon perturbations through 
Compton scattering. In turn, the baryon perturbations 
drive the behavior of the DDM perturbations, very ef- 
ficiently before DDM decoupling so that the DDM den- 
sity contrast 5 X on scales that enter the horizon during 
this epoch track the oscillations of the baryon-photon 
fluid before growing, and less efficiently after DDM de- 
coupling so that the baryons simply cause a drag on the 
growth of S x . In either case the matter power spectrum 
is suppressed relative to the standard case. The behav- 
ior of the CMB angular power spectrum can be simi- 
larly understood. Roughly speaking, the coupling of the 
DDM and baryons increases the effective baryon load- 
ing of the plasma at early times so that the CMB power 
spectra look similar to those from high-baryon models. 
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FIG. 12: Matter power spectra including baryon-DDM drag. 
The solid (red) curve is for is for (T> 2 + M 2 ) 1/2 = 1.4 x 
10" 15 e cm The short-dash (black) is for (V 2 + M 2 ) 1/2 = 
1.0 x 10 -16 e cm The long-dash (blue) curve is for is for (T> 2 + 
M 2 ) 1/2 = 5 x 10~ 15 e cm. These are all for a mass of 1 GeV. 
The curves are all for the standard concordance cosmological 
parameters, and the data points are from SDSS 42] 



m x and (D,A4) after marginalizing over other cosmo- 
logical parameters. We consider flat ACDM models and 
our chosen parameter space is the dark-matter density 
Q m h 2 , the baryon density flbh 2 , the Hubble parameter h 
in units of 100 km sec -1 Mpc -1 , the optical depth tcmb 
to the last-scattering surface, and the primordial spec- 
tral index n. We have employed the Markov chain Monte 
Carlo technique (see, e.g., Ref. |41j) to efficiently explore 
this parameter space, taking the most recent results from 
SDSS [13, WMAP |H, CBI 0, VSA gjj, and SNe la 
|46| as our data. Note that although DDM has no effect 
on observations of Type la supernovae, we include these 
data because the other parameters we allow to vary are 
constrained by these observations. We conclude using a 
relative-likelihood test that cosmological measurements 
lead to the bound shown in Fig. ^ The numerical calcu- 
lations confirm the qualitative behavior discussed above. 
Dipole moments as large as T> ~ 10~ 17 e cm, near the 
upper end of our allowed parameter space, are thus cos- 
mologically viable. 



VII. GAMMA RAYS 
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FIG. 13: CMB power spectra including DDM-baryon drag. 
The labeling of the curves is the same as in Fig. 1121 and the 
data points are those from WMAP |43 | 



This is of course an imperfect correspondence as modes 
of larger wavelength enter the horizon when the coupling 
is weaker, and so at later and later times the evolution 
of the photon perturbations becomes more and more like 
the standard-CDM case. But due to geometrical pro- 
jection effects modes of wavenumber k contribute to all 
I < kd,A where d,A is the angular-diameter distance to 
the last-scattering surface, and so the effects of DDM on 
small length scales can be noticed even on relatively large 
angular scales in the CMB. 

As the effect of DDM on the CMB may be partially 
degenerate with other cosmological parameters, we have 
explored a parameter space that allows us to constrain 



DDM particles in the Galactic halo can annihilate to 
two photons through the diagrams shown in Fig. |3J Since 
halo particles move with velocities v ~ 300 km sec" 1 <C c, 
the photons produced will be very nearly monoenergetic 
with energies equal to the DDM-particle rest mass. The 
intensity at Earth of such gamma rays is obtained by in- 
tegrating the emissivity, n 2 ,(a xx ^2'yv) , where n x is the 
DDM number density, along the given line of sight. The 
intensity is largest toward the Galactic center, where 
the dark-matter density is largest. In this direction, the 
gamma-ray intensity is then Q, 
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do = LOxl °" 



10 



r XX^2 7 ' 



10 -30 cm 3 sec" 
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(28) 

where / is a scaled integral of n 1 along a line of sight 
toward the Galactic center. The numerical coefficient is 
one-half that from Ref. [i[ since we have here particle- 
antiparticle annihilation rather than Majorana annihila- 
tion. Roughly speaking, I ~ 3 — 30 for cored-isothermal- 
sphere models of the Galactic halo, while / can extend 
up to ~ 300 for Navarro- Frenk- White profiles [42J; i.e., 
uncertainty in the dark-matter distribution in the inner 
Galaxy leads to an uncertainty of two orders of magni- 
tude in the predicted flux. We thus expect 
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3x 10" 11 ) 
17) cm " 2 



-M 4 



-1 „_— 1 



sec sr 



(29) 



To constrain dipole moments from non-observation of 
a gamma-ray line, we choose to use the most conserva- 
tive estimate, / ~ 3 for the dimensionless line integral. 
Moreover, we are not aware of any EGRET analysis that 
places limits in particular to a line flux. We thus obtain 



very conservative limits by using the binned continuum 
fluxes for the total gamma-ray flux listed in Table 1 of 
Ref. |4jj and noting that a line flux in that bin cannot 
exceed the measured continuum flux. The EGRET limits 
apply for masses 0.1 < (m x /GeV) < 10, and range from 



X>i7 < 180 for m x < GeV to V 17 < 100 for 



10 GeV, 



as shown in Fig. ^ 

Again, a few caveats are in order. First of all, our limit 
is quantitatively conservative — we chose the halo model 
that produces the lowest flux, and a detailed EGRET 
analysis would probably yield a line-flux limit lower than 
what we have assumed. On the other hand, the strong 
dependence ex T> 4 of the predicted flux on the dipole mo- 
ment guarantees that the upper limit to the acceptable 
dipole moment will not depend quite so strongly on these 
details. Second, if T>n > 5 in the mass range 100 MeV 
to 1 GeV, then the correct cosmological abundance most 
likely requires a particle-antiparticle asymmetry. If so, 
then the annihilation rate in the halo could be reduced 
far below the values we have obtained above. We con- 
clude by noting that with the increased sensitivity of the 
Gamma- Ray Large Area Space Telescope (GLAST), a 
detailed search for a line flux, and the possibility that 
the actual halo model provides a more generous annihi- 
lation rate, an observable GLAST signature may exist 
for masses 0.1 — 1 GeV and dipole moments as low as 
V 17 ~ 10. 



VIII. DISCUSSION 
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an m x < MeV particle has a dipole moment [D, M\ < 
10~ 22 e cm and no other interactions with ordinary mat- 
ter, then it might still be consistent with BBN. Of course, 
such a particle will, assuming standard freezeout, have a 
mass density many orders of magnitude larger than the 
dark-matter density. But suppose we were to surmise 
that the dark-matter density was fixed by some other 
mechanism. E.g., suppose the dipole was sufficiently 
weak that it never came into equilibrium. In this case, an 
additional constraint to the dark-matter dipoles can be 
obtained from energy-loss arguments applied to stars in 
globular clusters. Such arguments eliminate dipole mo- 
ments [T>, A4] <6x 10~ 23 e cm for masses m x < 5 keV 
|49J | . We have also considered constraints from astrophys- 
ical phenomena such as the stability of the Galactic disk, 
lifetime of compact objects, and annihilations in the so- 
lar neighborhood Q, and find that these constraints on 
the mass and interaction strengths are not competitive 
with those presented here. 

It would be of interest to attempt to embed this sce- 
nario in a consistent particle physics model. We might 
find links between baryonic and non-baryonic matter 
abundances, the dark matter electric dipole moment and 
the CP violation needed for baryogenesis, and the mag- 
netic moments of dark matter and baryons. However, 
such model building is beyond the scope of the present 
study. Our approach throughout has been entirely phe- 
nomenological, as we have been motivated by the desire 
to answer the question, "How dark is 'dark'?" 



In this paper we have considered the cosmology and 
phenomenology of dark-matter particles with a nonzero 
magnetic or electric dipole moment. We have found 
that information from precision tests of the standard 
model, direct dark-matter searches, gamma-ray experi- 
ments, and the CMB and large-scale structure restrict 
the dipole moment to be fDjA'f] < 3 x 10~ 16 e cm for 
masses m x < few GeV and [2?,^] < 10 _24 e cm for 
larger masses. Some of the allowed regions of param- 
eter space may soon be probed with GLAST and with 
future more sensitive direct-detection experiments. The 
electromagnetic interactions of these particles with nu- 
clei are coherent. Moreover, these particles cannot an- 
nihilate directly to neutrinos. Therefore, searches for 
energetic neutrinos from decays of the products of XX 
annihilation in the Sun or Earth are thus likely to pro- 
vide less sensitive probes than direct searches [48j . More- 
over, if there is a particle-antiparticle asymmetry, then 
the energetic-neutrino flux could be reduced without al- 
tering the direct-detection rate. 

We have restricted our attention to particles with 
masses m x > MeV, with the notion that lower-mass par- 
ticles will violate BBN limit, as discussed toward the end 
of Section |Hj We also consider masses m x > MeV, as 
particles of lower mass will almost certainly undergo rel- 
ativistic freezeout and thus lead to unacceptable dark- 
matter candidates. However, as also noted above that if 
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APPENDIX A: DERIVATION OF THE 
BARYON-DARK MATTER COLLISION TERM 

To determine how the cosmological perturbation equa- 
tions for baryons and dark matter are altered when we 
bestow the dark matter with a magnetic or electric dipole 
moment, we must formally evaluate the collision opera- 
tor of the general-relativistic Boltzmann equation in a 
given gauge |13L 13-H l35| for the dipole interaction of Eq. 
(Q. We have completed this calculation in detail, and 
find that the dipole interaction produces a drag force 
proportional to the relative velocity V — v x — Vb of the 
dark-matter and baryon fluids. As the relative velocity is 
gauge invariant in linear perturbation theory and all scat- 
terings are local processes, we may thus take a simpler, 
more physically transparent approach and just evaluate 
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this drag force using nonrelativistic statistical mechanics. 
It is this approach we now present. 

We wish to calculate the drag force per unit mass, or 
deceleration, due to collisions with protons to the dark- 
matter fluid as it passes through the baryon-photon fluid. 
Comoving scales A > Mpc enter the horizon when the 
cosmological temperature is T < 10 eV, when the DDM 
particles (which are restricted to m x > MeV) are non- 
relativistic. We may thus consider thermal velocity dis- 
tributions for nonrelativistic baryons and dark matter. 
Since the drag force can only depend on the dark-matter- 
baryon relative velocity, we take the baryon fluid to be at 
rest and the dark-matter fluid to have a velocity of magni- 
tude V in the x direction. Then, the proton phase-space 
distribution is 



fp(v P ) = 



(27r) 3 / 2 i;3 



-</2vl 



(Al) 



where v p = (kTp/rrip) 1 / 2 is the proton velocity dispersion 
and n p the proton number density, and 



fx@x) 



(2^) 3 / 2 w3 



exp 



(v x - Vx) 



21 



2v x 



(A2) 



is the dark-matter phase-space distribution, with v x = 
{kT/m^ 1 / 2 . Recall also that we expect V <C v p , as 
discussed above. 

The drag force per unit mass is obtained by integrating 
the momentum transfer per collision over all collisions 
between protons and dark-matter particles. From the 
symmetry of the problem, the deceleration of the dark- 
mattcr fluid will be in the x direction, and it will have a 
magnitude, 



— / d 3 V x f x {v x ) J d 3 Vpfp(v p )\v p - v x 



dn^(v xxf 



J X xl ' 



(A3) 



Here Vt = (9, <fi) is the scattering angle in the center-of- 
mass frame, and v xx f — v XX i is the difference between 
the final and initial x component of the dark-matter- 
particle velocity; the difference is the same in the center- 
of-mass and laboratory frames. The differential cross sec- 
tion da/dfl is that given in Eq. (J5J. 

Consider an individual scattering event. Let a be the 
angle that v p — v x makes with the x direction; this is then 
the angle that v x makes with the x axis in the center- 
of-mass frame, and the magnitude of the initial and fi- 
nal dark-matter velocities in the center-of-mass frame 



is u™ 1 = m p v/(m p + m x ), where v = \v p 



v x \ is the 



relative velocity. The initial x component of the dark- 
matter velocity in the center-of-mass frame is then v XX i = 
v x m cosa. The scattering angles 9 and </> are then the po- 
lar and azimuthal angles that the scattered dark-matter 
velocity make with the initial velocity in the center-of- 
mass frame. By rotating this coordinate system by an 
angle a about the z axis to align it with the laboratory 
x axis, we find v xx f — 
Thus, 



1 (cos a cos 9 — sin a sin 9 sin < 



/ dn ~d~&( v xxJ ~ v x*i) 



m p Z 2 e 2 V 2 cos a 
27r(m p + m x )v 



(A4) 



Completing the integral in Eq. I|A3[) in the limit V <C 
(v p ,v x ), we find 



2 (ZeV) 2 m p n p V 



37r(2 7 r) 1 /2 mp + 



! ^ —I— 1 ! ^ 



(A5) 



v ^ " b x Jvp]+v x 



Taking into account the definition (v p ) = x /8T p /(7rm p ), 
this drag force leads to the drag- force term in Eqs. {T7J|, 
(|20|l . and (|21() . when including the simple corrections for 
a mass fraction Y of helium. 
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